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Abstract 



(N 

CN \ We consider the dynamics of a relativistic Dirac particle constrained to move 

in the interior of a twisted tube by confining boundary conditions, in the 
J> ■ approximation that the curvature of the tube is small and slowly varying. In 

contrast with the nonrelativistic theory which predicts that a particle's spin 
does not change as the particle propagates along the tube, we find that the 
angular momentum eigenstates of a relativistic spin-^ particle may behave 
£^ ' nontrivially. For example, a particle with its angular momentum initially 

polarized in the direction of propagation may acquire a nonzero component 
of angular momentum in the opposite direction on turning through 2ir radians. 
Also, the usual nonrelativistic effective potential acquires an additional factor 
<D ' in the relativistic theory. 

43 



I. INTRODUCTION 

Twisted tubes, in which particles are confined by an external potential, have proven to be 
excellent settings for the study of quantum mechanics in curved geometries. They have been 
studied extensively in the context of nonrelativistic quantum mechanics, with particular fo- 
cus on the case of small and slowly varying curvature Jl|. One important example is the 
limit of a very narrow tube, where the "confining potential" constrains a particle to move on 
a curve, just as a physical bead might be attached to a wire. Subsequent authors have con- 
sidered the natural extension to a particle confined to a small neighborhood of an arbitrary 
Riemannian manifold embedded in a higher-dimensional Euclidean space and have applied 
their formalism to the energy spectra of molecules ||. For a recent survey of the dynamics 
of a particle in various confining potentials in a loosely interpreted "tube" see Ref. ||. How- 
ever, analysis has been largely limited to nonrelativistic quantum mechanics, in which spin 
does not affect the dynamics; the relativistic theory has been thoroughly investigated only 
for confinement and quantization schemes other than the confining potential approach [|J. 
Here we present for the first time a treatment of a relativistic Dirac particle in the confining 
potential approach using confining boundary conditions. 
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The phenomena associated with the propagation of a relativistic fermion within a tube 
contrast with the simpler and well-known case of light propagating in an optical fiber, 
studied in the mid-1980's by Chiao and Wu 0. They assumed the fiber to have constant 
circular cross-section and small and smoothly varying curvature. In this approximation, 
they argued that the two states of circular polarization (helicity ±1) propagate along the 
fiber without mixing with one another and accumulate opposite Berry phases proportional 
to the integral of the torsion along the fiber. It is instructive to re-examine their argument 
that the helicity eigenstates do not mix g]. First they use the adiabatic theorem to argue 
that modes with different transverse waveforms do not mix as they propagate through a 
smooth fiber. The lowest energy level is doubly degenerate, with total angular momentum 
component along the axis of the fiber J ■ t = ±1. The longitudinal mode, J • t = 0, is 
excluded by gauge invariance. Following Chiao and Wu we refer to these modes somewhat 
loosely as photon states of helicity ±1. In the limit of a thin fiber, the only surviving terms 
in the Hamiltonian affecting helicity are of the form k ■ S, where the "vector curvature" , k, 
is a vector in the direction of the principal normal to the fiber and of length k, 1/k is the 
instantaneous radius of curvature of the fiber, and S is the photon's spin. Since k ■ S is a 
vector operator, the Wigner-Eckart theorem allows a change in helicity only by Ah = ±1. 
Ah = ±2 is forbidden, so the helicity states propagate without mixing, even though they 
are degenerate - a consequence of gauge invariance. 

The point of departure for our paper is the observation that Chiao and Wu's argument 
does not apply to a Dirac fermion, which has doubly degenerate modes with helicity J -t = 
±(n+ ~). Once again we refer to these modes somewhat loosely as fermion states of helicity 

±(n + |). In particular, the ground state modes have n = 0, so J • t = ±|. The Wigner- 
Eckart theorem does not forbid transitions between these states because the eigenvalues of 
J ■ t differ by 1. We cannot exclude the possibility that the degenerate angular momentum 
states will mix as they propagate on the basis of symmetry. Therefore the Dirac equation 
demands more careful analysis. 

As in Ref. |l[ , we consider a particle confined to the interior of a tube with circular cross- 
section. We work in the rest frame of the tube. The tube may be constructed by transporting 
a circular disc such that the center of the disc traces a smooth curve X(s) and the tangent to 
the curve is normal to the disc for all values of the arc-length parameter s. The region swept 
out by the disc is the tube. We model the potential that confines the particle to the tube by 
means of a boundary condition. The boundary condition must guarantee that the normal 
component of the probability current j = ip^aip vanishes on the surface of the tube. In the 
nonrelativistic problem one can simply impose Dirichlet boundary conditions. However, in 
the Dirac theory, the requirement that all components of the wavefunction vanish on the 
boundary is too stringent, for then in general the Dirac wavefunction will vanish everywhere. 
A less severe condition that will still confine probability is 

— ij ■ nip — ip (1) 

where n is the unit normal (in the outward sense) to the boundary of the tube 0. That 
eq. (P guarantees probability confinement follows from the commutation relations for the 
7-matrices and the antihermiticity of 7. Eq. ([TJ) reduces to the usual Dirichlet boundary 
condition on the "large" components of the Dirac spinor in the nonrelativistic limit. 

It is a general feature of the Dirac equation that spin and orbital angular momentum 
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are not separately conserved. It is not surprising, therefore, that the boundary condition, 
eq. (|l|) does not commute with the spin operator. The result is the nontrivial transport 
of angular momentum states that we will calculate. We could have chosen other boundary 
conditions, for example the chiral generalization of eq. ([[]), 



parameterized by the chiral angle 6, with a slightly different phenomenology. We will explore 
some of the alternatives in the final section of this paper. 

In Section II we establish coordinates in the interior of the tube that will be convenient 
for our purposes. Section III is a construction of the dynamical equations of our system in 
the parallel transport frame. We write the dynamical equations as those for propagation in 
a straight tube with an interaction. In Section IV we calculate the necessary basis states 
for a Dirac particle in a straight tube. In Section V we specialize to the approximation that 
the curvature is small and slowly varying. This approximation appears to capture most of 
the interesting physics. We then expand the Hamiltonian to lowest nontrivial order in small 
parameters. Using the resulting Hamiltonian, we compute in Section VI the relevant dynam- 
ical equations. Section VII contains our conclusions and also briefly discusses possibilities 
for boundary conditions other than eq. ([I]). 



We construct a coordinate system for the tube described in Section I by erecting a right- 
handed orthonormal frame at each point on the curve. We parameterize the curve as X(s) 
where s is the arc-length parameter. Since the tangent to the curve is always normal to the 
cross-sectional disc that sweeps out the tube, a natural choice of frame is {ei, £2, £3} where 
t\ and t2 are orthogonal unit vectors in the plane of the disc and £3 is the unit tangent to 
the disc in the direction of increasing s. For clarity we adopt a summation convention where 
Greek indices, /1, u, etc., are summed over 1, 2, and 3, whereas Latin indices, i, j, etc., are 
summed only over the "transverse" directions 1 and 2. To completely define the frame, we 
must specify the s-dependence of e\ and e 2 . We require that the frame rotates smoothly as 
the disc sweeps out the tube. That is, 



where uj(s) is a smoothly varying instantaneous "angular velocity" associated with the ro- 
tation of the frame. 

Since the curve uniquely determines e 3 (s), it determines u(s) up to the addition of a 
constant multiple of e^(s), corresponding to an additional rotation of the transverse axes 
about the tangent. For simplicity, we adopt the "parallel transport" frame, in which u has 
no component along the tangent. There is a simple visual model of the parallel transport 
frame. Imagine a rigid wire with a disc attached to it by a frictionless sleeve that keeps the 
disc normal to the sleeve. Without friction all the torques act in the plane of the disc so 
the angular momentum about the tangent direction is conserved. If the disc has no angular 
velocity about £3 initially, it will not develop any. In this case, vectors drawn on the disc 



— i'y ■ hifj = e* 75 ip 



(2) 



II. ELEMENTARY FORMALISM 




(3) 
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are said to execute parallel transport along the wire. A simple calculation shows that the 
instantaneous angular velocity of this frame is 

- = e 3 x — . (4) 

A new coordinate system for points inside the tube may now be specified in terms of the 
parallel transport frame by 

In this coordinate system, the spatial metric tensor is diagonal 

/I \ 



1 

V (i-s-SV 



(6) 



where the vector curvature k of X(s) is defined as 

de 3 



h h — — = k 1 ^! + /t 2 e 2 . (7) 
as 



Comparison with the definition of uj shows 



uj = e 3 x k or 

u 2 = k and to 1 = —k 2 . (8) 



The diagonal metric tensor yields a simple gradient operator: 



~ 9 . d A 1 9 

ei^-r + ea^ + eg — . (9) 



The volume element for integration in these coordinates is y det gj^dsd 2 ^ = (1 — k ■ ^)dsd 2 ^. 

Since the tube intersects itself if k ■ £ equals one, we require that kR < 1 where R is the 
radius of the disc. 

An important subtlety is that while the choice of the parallel transport frame is always 
valid locally, it is not always possible to adopt it globally. In particular, if the curve is closed, 
parallel transport of the frame through one circuit of the curve need not return the frame 
to its original orientation. This geometric feature manifests itself in classical physics as the 
well-known Hannay's angle [f§|l]. For the purposes of this paper, we assume that our tube 
continues indefinitely without intersecting itself. 

III. DIRAC EQUATION IN THE PARALLEL TRANSPORT FRAME 

In free space, with units such that h = c = 1, the Hamiltonian in the Dirac equation 
m> = is 

U = (-ia ■ V + l3m + V) . (10) 
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We define a and (3 as usual by a 1 = 7 7' and (3 = 7 ; in the chiral basis that we employ. 



In the coordinate system defined in Section II, the Dirac Hamiltonian takes the form 

d i 

H = -ia'e j -Q-j--a-e 3 fo+Prn (12) 

where A = 1 — k ■ £. We normalize our wavefunctions \I> so that 

/" cPx&V = J dsd 2 £A^^ = 1 . (13) 

The factor of A in the integration measure complicates calculations and questions of her- 
miticity. Therefore we define new wavefunctions i/j by 

%P = aH (14) 

which have the normalization J dsd 2 ^^ip = 1. The original Dirac equation can be re- 
cast in terms of a Hamiltonian operator H, defined by H = AzTiA^^, which acts on the 
wavefunctions ip, 

.-> - 9 i _ _ i _ A rd 1 d{K-U\ _ 

H = —ta-ej—— -a • k — —a ■ e% — — I + pm . (15) 

3 d& 2A A ^ds 2A ds > H v ' 

This Hamiltonian is difficult to calculate with because the matrices a-e^ are s-dependent. 
That is, the momentum operators are expressed with respect to the rotating coordinate 
frame while the spin operators are expressed with respect to some fixed rectilinear frame. 
To remedy this problem, we change the spin basis. We define a unitary transformation Q 
by 

Q(s)a- ^(3)^(3) =« M (16) 

for all s. Here we have made the s dependence of the unit vectors e M explicit. 

Differentiating the above condition with respect to s and using eq. (|3|) for de^/ds, we 
find 

^ = -S-c^ (17) 
ds 2 K ' 



where S* are the 4x4 Dirac spin matrices. Therefore, 



Q(s) = V exp 



(18) 



where V denotes path-ordering along the tube. 

To take advantage of this transformation, we once again define new wavefunctions, 



x(s,() = Q(s)ij(s,i) (19) 
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and a new Dirac Hamiltonian, H, 

H = QHQ^ . (20) 

The form of eq. fllBp indicates that the transformation Q(s) just represents a change of spin 
basis to the spinors representing spin up and down in the e 3 (s) direction. Although this 
transformation replaces the a-matrices as intended, it generates a "gauge term" in H, 

Simplifying a 3 £ - a; yields —id ■ k. Thus, this gauge term cancels the second term in eq. fllBD 
giving the following relatively simple Hamiltonian for \\ 

H = _ !Q 3_ + l3m _ la ,(_i_ + _ A_M) . (22) 

The transformations we have performed reduce our problem to the study of this perturbed 
Dirac Hamiltonian in a right circular cylinder, (£ 1 ) 2 + (£ 2 ) 2 < R 2 , —oo<s<oo, with respect 
to a fixed rectilinear frame. The price of this simplification is the variety of curvature 
dependent terms in eq. (p2|). 



IV. DECOMPOSITION INTO TRANSVERSE EIGENSTATES 

We look for eigenstates of the Hamiltonian in eq. (p2|) of the form 

x(s,£) = Y, c {p}( s Mp}(£) (23) 

M 

where the {?7{ p }} are eigenstates of the transverse Hamiltonian, 

H ± = -ta^ + f3m (24) 
satisfying the boundary condition 

-il-rr] = rj\ I=R . (25) 

Here we have introduced plane polar coordinates, r = |£| and ip = tan -1 ^ 2 /^ 1 . The labels 
{p} denote the quantum numbers necessary to specify the eigenstates of H± completely. 

H± is just the free Dirac Hamiltonian in one lower dimension. 5*3 and L3, the tangen- 
tial components of the spin and orbital angular momentum respectively, do not separately 
commute with H± but J3 = L3 + 5*3 does. J3 also commutes with the operator 7 • r in the 
boundary condition. In addition, the operator 7°S , 3 commutes with H±, J 3 , and 7 ■ f. We 
define the "transverse energy" as the eigenvalue A of H±. An eigenstate of H± with eigen- 
value A may thus be chosen to be an eigenstate of the helicity operator J3 with eigenvalue v 
{y — n — 1/2 for some integer n) and of the operator 7°S , 3 with eigenvalue e (e = ±|)- For 
simplicity, we use the notation {v^, A} when specifying these quantum numbers. 
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H±_ has the usual Dirac charge conjugation symmetry, so each eigenstate with eigenvalue 
+ |A| is paired with another with eigenvalue — |A|. Since a 3 anticommutes with H± and 7°S , 3 
and commutes with J 3 and j-r, a ir q u ± x = r\ v T-\. The transformation that sends J3 to — J3 
commutes with H± and 7 • f and anticommutes with 7°S , 3. Thus, the states with z/ + have 
the same eigenvalue A as states with —v~. For each value of v and e, there are towers of 
radial excitations, one for positive A and one for negative A, obtained by adding more nodes 
to the radial wavefunction. A summary of the eigenstates of H± for small v and A is given 
in Fig. 0. 



XR 
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FIG. 1. Summary of the low-lying transverse energy eigenstates in the extreme relativistic limit, 
m — > 0. Only positive A are shown. Each level is twofold degenerate corresponding to v + <-> — v~ . 
Radial excitations are displayed with dashed lines to correlate with Fig. 0. 



With these general considerations in mind, we construct the explicit representation of 
the wavefunction rj u ± t \ for positive A: 



-iJ n (qr)e inip 
±( ±1 J n -i(qr)e i{n - 1 ^ 
±i.J n (qr)e inip 



(26) 



where n is an integer with v — n— J n is a Bessel function of the first kind, q = +VA 2 — m 2 , 
^ _ m+^/q +m _^ ^ satisfies the following equation obtained from the boundary condition: 



J n {qR) = ±C ±l Jn-i{qR) ■ 
The normalization N is fixed up to a phase by the condition 

d 2 ^rfr] = 1 . 



(27) 



(28) 



The unusual form of the spinors in the nonrelativistic limit is due to our choice of the 
chiral representation of the Dirac matrices (75 = diag(— 1, — 1, 1, 1)). Transforming the 
eigenstates to the Dirac basis (70 = diag(l, 1, — 1, — 1)), we find that only one of the two 
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upper components and one of the two lower components of each wavefunction are nonzero. 
In the nonrelativistic limit, for the states with positive A, the lower components in the Dirac 
basis tend to zero. These are eigenstates of 5*3 and L 3 , in agreement with the fact that 
7°S , 3 is just 5*3 on the upper components in the Dirac basis. We use the chiral basis because 
the wavefunctions can be expressed more compactly away from the nonrelativistic limit and 
because computation of matrix elements is simpler. 

For a straight tube, the energy eigenstates with helicity v and longitudinal momentum 
k are obtained from the ansatz, eq. dffi): 



k sgn(A) 



V\ 2 + k 2 + \X 



(29) 



In a twisting tube, the energy eigenstates could be superpositions of several {^±,a} with 
nontrivial {c u ± : \(s)} determined by eq. (|22|). 

In later sections, we will study the propagation and mixing of the ground states of 
the transverse Hamiltonian, i.e., those with the smallest possible value (Ao) for |A|. For 
concreteness, we present the explicit form of the wavefunctions for this fourfold degenerate 
set of states: 



N 



V2 = ?L.i- 

2 



,A 



N 



( (Mqr) \ 
-U x {qr)e l * 

(Mqr) 
V U x {qr)e^ J 



( -U 1 (qr)e- i ^\ 
iJ x (qr)e~ llp 



V3 = Vl 



-Ao 



N 



( (Mqr) \ 

U x {qr)e l * 

-CJo(q r ) 

\iJ x {qr)e ilp J 



r)4 = -rj_i + 

2 ' 



-Ao 



N 



(U 1 {qr)e- iv \ 

C,Mqr) 

V -(Mqr) J 



(30) 



where we have moved some of the factors in the representation in eq. ( pf) into the normal- 
ization for later convenience and used J_ n = (—l) n J n . The eigenvalue condition 



MqR) = ^Ji(qR) 

determines Ao, and the normalization condition eq. ( 

1 _ _o/(C 2 + l) 2 2 



(31) 



N 2 



2nR 2 (- 



c 2 



CqR 



then fixes N: 



Jl{qR) . 



(32) 



In the nonrelativistic limit (m ^> q), eq. ( pT|) reduces to Jo(qR) = 0, a familiar result, and 
in the ultrarelativistic limit (m — > 0) it becomes Jo{qR) = J\{qR). 



V. ADIABATIC CONDITION 

In addition to assuming; that the vector curvature, smooth function of s, we 

assume that it is nonzero only in a finite region. Thus the tube consists of two infinitely 
long straight regions joined by a curved region of finite length. Our goal is to propagate the 
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unperturbed solutions to the Dirac equation, labelled by their longitudinal momentum, k, 
and their transverse quantum numbers, {z/^A}, through the region of nonzero curvature. 
We expect both scattering states (k real) and bound states (k imaginary). By the "transverse 
quantum numbers" of a tube eigenstate, we mean the quantum numbers associated with the 
first term of its representation in the form given by eq. ( p9"D . 

We restrict our analysis to the case of adiabatic propagation in a mildly curved tube 
that dominates the work in this field. This amounts to the conditions 

\kR\ < 1 

< 1 . (33) 

In this approximation to the nonrelativistic problem, states with different transverse wave- 
forms do not mix ||. Because of the analogy between time development with a slowly 
varying V(t) and s-development with a slowly varying curvature, this is known as the "adi- 
abatic approximation" . The generalization of the adiabatic approximation to the Dirac case 
is that only states with the same values of |A| can mix when propagating in a smooth enough 
tube. For example, the four states with \v\ — ^ given explicitly in eq. (|30|) can mix adiabat- 
ically. Since this subspace includes states of different J • £3, interesting angular momentum 
transport is still possible in the adiabatic limit. 

For simplicity, we consider the case of small longitudinal momentum. We assume that 
the square of the longitudinal momentum of a state under consideration is much smaller 
than the difference between the square of its transverse energy and that of other states with 
differing transverse energy. For a particular state with a set of transverse quantum numbers 
p , this amounts to the condition 

\kR\ <^d p (34) 

where d p = min p / ^J\(q p iR) 2 — (q p R) 2 | and p' runs over all sets of transverse quantum numbers 
such that q p i 7^ q p . 

For a physical tube, these conditions are easily satisfied. For a fixed mass m and a 
given curve with smoothly varying vector curvature k(s) which is nonvanishing only in a 
finite region, we can always choose a nonzero R small enough to satisfy each of the above 
conditions. There are, however, some subtleties associated with this condition. For a fixed 
tube radius R, making m arbitrarily large makes d p arbitrarily small for most p. The 
additional degeneracies that develop as mR — > 00 originate in the decoupling of spin and 
orbital angular momentum. Thus, pairs of states with L 3 = 1 and S3 = ±|, which are 
split by relativistic effects, become degenerate as mR — > 00 as shown in Fig. (|2]). The cases 
p — {- , A} or { — I , A} are an exception. These states, which have L 3 = and S3 = ±~ 
are degenerate for all mR, and no new degeneracies arise as mR — > 00. Again, see Fig. (0). 

Of course, this is not problematic because we know exactly how spin and orbital angular 
momentum eigenstates behave in the nonrelativistic limit: the Dirac problem reduces to the 
scalar problem solved, for example, in Ref. [|l|, accompanied by a dynamically trivial spin 
label. 

Each of the conditions in eqs. (|3^) and (0) has a small parameter associated with it. 
We let S denote the largest of these parameters. To simplify the analysis of the relative 



da 3 



ds 
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2 4 6 8 10 

mR 

FIG. 2. Dependence of qR on mR for the five lowest levels in Fig. [I] 

importance of terms in the Hamiltonian, we introduce scaled coordinates u? = £ J /R and 
multiply the Dirac equation by R to obtain H'x = E'x where H' = RH and E' = RE are 
the scaled Hamiltonian and energy respectively. From eq. (|22|), we obtain 



TT , .id n _ . ,_c? (RK-u^d R d(RK-u)\ 

So far, this equation is exact. The first two terms are zeroth order in 8, the third term 
is first order in 5, and the last two terms are second order in 5. To second order in 5, we 
can replace A -1 by 1: 

H' = -ia j ^L + f3mR - ta 3 ((l + Rk ■ u) R^- + g ^gLjh + (8 3 ) . (36) 
ou J v as 2 as ' 

For later use, we also compute H' 2 to second order in 5: 

H» = -Vl + - - • -| - f | " -) + 0(^ 3 ) • (37) 

Here E is the Dirac spin matrix, S = ( ^ ^ 

V a. 

Having systematically expanded the Hamiltonian to second order in small quantities, we 
may now return to unsealed coordinates. The adiabatic condition implies that the effective 
potential mixes two transverse eigenstates rj p and rj p > to second order in 5 only if |A P | = |A p /|, 
or equivalently if q p = q p >. Thus, to second order in S, the mixing amplitudes outside of the 
four dimensional subspaces of fixed q are negligible. 

As is often the case with the Dirac equation, it is easier to work with the square of 
the Dirac operator instead of the first order form. Since the curvature vanishes for s < 
s for some s , the solution x( s ^0 is given by the straight cylinder problem for s < s . 
This determines both x( s o>£) an d dx(s, £,)/ds\ s=SQ . To solve the Dirac equation to the 
desired order, we project the equation onto the four- dimensional subspace containing x( s o> I) 
to obtain four coupled first order equations. Equivalently, any solution of the projected 
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equation must also satisfy the projection of the equation H 2 x = E 2 x to second order in 5. 
The initial conditions on x an d dx/ds at s = So determine a unique solution of the second 
order equation, which must also be the unique solution of the first order equation. From 
now on, we study the second order equation, eq. (|37D . We complete the square and divide 
by R 2 , to obtain 

H 2 = -V? + m 2 - - -t ■ uf - ^ + O^/R 2 ) . (38) 

Notice that the terms of order 5 3 /R 2 , can be neglected in the adiabatic approximation. H 2 
is simpler to study than H because it only mixes states with the same A. Thus the mixing 
problem has been reduced from a4x4toa2x2 problem. For simplicity, we refer to the 
two dimensional space of states (with J ■ t = ±v) that can mix under H 2 as the "adiabatic 
subspace" . [This simplification is analogous to the Foldy-Wouthuysen transformation that 
eliminates the lower components of the Dirac spinor in the nonrelativistic limit.] This 
remarkably simple result is interpreted in the following section. 

VI. PROPAGATING MODES AND SPIN EVOLUTION 

This Hamiltonian (squared) of eq. (|37| ) appears to be the same one that one would obtain 
in a nonrelativistic theory. The curvature induced "potential", — k 2 /4, is well known. The 
"gauge term", — |E ■ uj, looks like it undoes the rotation, fl(s), that put us into the moving 
frame. If it did, a stationary experimenter outside the tube would observe a spin vector that 
always points in the same direction as the particle propagates in s. However, this description 
of propagation by this Hamiltonian must not be correct: Because S • u fails to commute 
with the boundary condition eq. (P, it cannot in general be diagonalized with respect to 
the transverse basis states identified in Section IV. 

To obtain the correct physical description, we must remember that H 2 is to be projected 
onto the adiabatic subspace containing two states of opposite helicity {u = ±|). The 
externally prescribed vector curvature, k, explicitly breaks rotation invariance and allows 
states of different helicity to mix. However, the symmetry violating terms in H 2 are at most 
vector operators (note (S • uj) 2 = k 2 ), and by the Wigner-Eckart theorem can only induce 
transitions between states with helicity differing by unity. Thus mixing of different helicity 
states can only occur if \u\ = |. For all the states with \u\ > |, helicity, and in fact all 
transverse quantum numbers, are conserved. This is exactly analogous to the propagation 
of photons discussed by Chiao and Wu ||. 

States with \u\ = \ mix as they propagate. The crucial feature that distinguishes the 
relativistic and nonrelativistic situations is that the states with v = ±~ are not eigenstates 
o/S 3 . Instead they are eigenstates of J 3 and 7°S , 3 . Since the ground states do not acquire 
any additional degeneracies in the nonrelativistic limit, we can apply our formalism to them 
in this limit. The ground states become L3 = and S3 = ±| states. The gauge term in 
eq. ( PS] ) exactly undoes the rotation Q(s) and leads to trivial spin transport along the tube. 
However, in the relativistic case, the matrix elements of S between states with v — ±| 
are reduced by a "depolarizing factor" which reflects the mixing of spin and orbital angular 
momentum. 
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We study the spin transport in the interesting case of v = ±|. For concreteness, we 
calculate the matrix elements of the operator E • u in the lowest adiabatic subspace, {rjj}, 
enumerated in eq. (|30D 

{j\E-<2\k) = J 'd^tE-u;^. (39) 

For positive A, only the (12) and (21) matrix elements are nonzero, and for negative A, only 
(34) and (43) matrix elements are nonzero, 

J d 2 £ rj{ E • u r) 2 = / d 2 ^i]l E ■ u rj 4 

= 2N 2 (u l - iuj 2 )( 2 [ R rdr dtp J (qr) 2 
Jo Jo 

= 2nN 2 R 2 (u 1 - iu 2 )(( 2 + 1) Ji{qR) 2 

= (uj 1 - iu 2 )D (40) 
where we have defined the "depolarizing factor" D as 

d = ^ + a (4i) 

A similar expression holds for each excited subspace with \ v\ = ~. For a subspace containing 
a state in the tower over 771+ Aq , one just changes q accordingly in the expression for D. For 

a subspace containing a state 771- ., one needs to change q and replace ( by — C -1 , in the 

2 > A 

expression for D. The general expression is 

C ±2 (r ±2 + 1) 

D(q, ±) = ^ l ^ + ' (42) 

where the sign ± is the product of the signs of the transverse quantum numbers and is the 
same for the 4 basis states in each adiabatic subspace. An important caveat is that this 
formula fails to hold in the nonrelativistic limit for the lower choice of sign. The states with 
v — ~ and those with v — | both reduce to L — 1 and have the same transverse energy 
in the nonrelativistic limit; the additional degeneracies signal that these subspaces do not 
evolve adiabatically. However, as remarked earlier, we know that in this limit D — 1 and 
a straightforward calculation extending to a 16-dimensional subspace of states verifies this 
fact. 

The factor (uj 1 — iu 2 ) in eq. fl40|) shows the parallel between the matrix elements of 
E ■ uj in the {rjj} basis and its matrix in the basis of eigenstates of E3. The constant 
of proportionality, D, measures the rate at which the spin vector rotates in the parallel 
transport frame as the wave propagates. In the nonrelativistic limit, £ — > 00 so D — > 1, 
and the spin completely decouples from the evolution of the tube. If D were zero for some 
choice of parameters, the spin would remain fixed in the parallel transport frame as for the 
photon. In the ultrarelativistic limit, ( = 1 and qR ~ 1.434, the lowest root of the equation 
Jo(qR) = J\(qR). Substitution gives D ~ 0.767, so that even in the extreme relativistic 
limit helicity is not conserved. 
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We now project the second-order Dirac equation, H 2 \ = E 2 x, onto the lowest multiplet, 
{rjj}. Writing \ — 2w=i c j( s )Vj(Qi the 4x4 matrix equation for the four component vector 
C(s) = (ci,c 2 ,c 3 ,c 4 ) is 



/ d i ^ - a 2 9 k 2 D 2 



+ mf 



C{s) = E 2 C{s) 



(43) 



where T is the vector of (4 x 4) Pauli matrices, T 



that generates the rotations 



of the J3 = ±h eigenstates. Just as the Pauli matrices, Sj, rotate S3 



r 

T. 

,,. ,.,„.«,,,., , ■_■ — 1/2 states into 
S3 = —1/2 states, the matrices 2} rotate the J 3 = 1/2 states into J 3 = —1/2 states. 

Eq. fl4"3"| ) is our fundamental result. Relativistic effects appear in two places. First, as 
pointed out earlier, the rotation back to the laboratory fixed frame is only complete in 
the nonrelativistic limit, when D = 1. As q increases relative to m, D decreases. The 
ideal case of D = 0, in which the helicity is conserved, is not reached even in the extreme 
relativistic limit. The following section discusses opportunities presented by other choices 
of boundary conditions. Second, the depolarizing factor appears in the effective scalar 
potential, — k 2 D 2 /A. If we apply the obvious unitary transformation Cd(s) = Q^^C^s) 
(see eq. (Il8|) ), where 



D / ds'T-uj(s') 



£Id(s) = Pexp 



then the Dirac equation for C^(s) reduces to the scalar case 



(44) 



d 2 



k 2 D 2 



ds 2 



C D = k 2 C D 



(45) 



This equation can be solved numerically for a given tube, with the usual phenomenology of 
bound and scattering states. 



VII. DISCUSSION 

It is helpful to interpret our results first in the nonrelativistic limit. We expect to find no 
coupling to the spin because the nonrelativistic boundary conditions do not involve spin at 
all and the exact Hamiltonian commutes with S. It is reassuring to see that our formalism 
recovers this result: in the nonrelativistic limit, D — > 1, so that in the parallel transport 
frame, the spin appears to rotate with angular frequency — uj, where the minus sign indicates 
that the sense of the rotation is opposite that of the coordinate frame. In other words, the 
spin vector maintains its orientation in the three-dimensional Euclidean space in which the 
tube is embedded (which corresponds to the coordinate frame of some three-dimensional 
Euclidean experimenter). Away from the nonrelativistic limit, D 7^ 1, so the variation of 
spin is no longer exactly opposite to the rotation of the parallel transport frame. 

Now if D were to go to zero, the vector J would not rotate relative to the axis of the 
tube. This would parallel the phenomenon described by Chiao and Wu. However, D is never 
zero in the lowest energy subspace. Even for m = 0, D m 0.767. We obtain a complete 
description of D as a function of mR by solving for the lowest root of the eigenvalue condition 
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eq. ( p7|) as a function of m and substituting in the expression eq. ([IT]) for D. A plot of D 
versus mR appears in Fig. (|3|). The states with p = {| + , A} do have D — > 1 as m — > oo, 
as we expect. However, the states for which p = {\ , A} have D — > in the nonrelativistic 
limit. This "unphysical" result appears as an artifact of the adiabatic approximation, which 
breaks down for these states in the nonrelativistic limit, as explained in Section VI. 



0.6 

D 

0.4 



0.2 



2 4 6 8 10 

mR 

FIG. 3. Depolarizing factor D as a function of mR for the ^ + and ^ towers of states in Fig. |IJ. 

The effect of the factor D on the vector J has a simple geometric interpretation. As an 
example, consider a tube whose defining curve lies in a plane and such that it has two straight 
sections, \s\ > s , joined by a section with constant curvature k and center of curvature 
P. Q In addit ion, construct a complementary tube parameterized by the same arc-length 
parameter s satisfying the following conditions: it is parallel to the first for s < — sq\ it has a 
section of constant curvature Dk with center of curvature P for \s\ < sq; and it has another 
straight section for s > Sq. The smaller curvature implies a larger radius of curvature. The 
effect of the depolarizing factor can be summarized by the following construction: as the 
Dirac particle traverses the curved section of the tube, the vector J remains fixed relative to 
the tangent of the complementary curve at the corresponding value of s. For m = 0, if the 
tube curves through A0 = 2tt, then J falls behind the tangent by approximately 1.534 it. 

In our formalism, we have used the boundary condition eq. (|TJ) to simulate a hard-wall 
potential barrier. Many of the qualitative features of our results follow from the fact that spin 
is not a good quantum number for the Dirac equation subject to such a boundary condition. 
It is natural to ask how sensitive these results are to the specific choice of confining boundary 
condition. An obvious generalization given in eq. (H) is characterized by a chiral angle 9. We 
have chosen 8 = throughout; let us now take 9^0. If we set f = e^^ip, we can map 
this version back to the original boundary condition at the cost of introducing a "chiral" 
mass term into the Dirac equation: 



Actually, we require the transition from zero to nonzero curvature to occur smoothly over a long 
enough interval that the adiabatic approximation is applicable. 
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(a ■ p + (3me- il5e ) xp 1 = Eip' (46) 
-iy ■ rip' = ip' . (47) 

Repeating the calculations in the earlier sections shows that the sole effect of the chiral 
angle is to replace m by mcosO in the expression for the parameter ( and thus implicitly 
in the boundary conditions and depolarizing factor. Thus the analysis above, with an 
effective mass mcosO, applies to this family of problems. In fact, for cos# > 0, the plot 
in Fig. (Q) describes the depolarizing factor as a function of mR cos 6. When mRcosO 
becomes negative new phenomena arise. In particular, the boundary condition eq. (D) acts 
like an effective attractive potential at r = R. As a result for sufficiently negative mR cos 9 
a qualitatively new solution appears in each partial wave, which is "bound to" r = R and 
decays exponentially as r decreases. The detailed study of these states is beyond the scope 
of this paper. 

In addition to these boundary conditions, there are other families of boundary condi- 
tions that confine the probability current but which we have not yet considered in detail. 
Specifically one can generalize eq. (H) by the condition 

— iy ■ nip = e r tp (48) 

where T is any antihermitian Dirac matrix that anticommutes with y ■ n. V need not be 
constant, as the demonstration of probability confinement involves only algebraic manipu- 
lations. Interesting alternatives include T\ = y 5 6(s), T 2 = 7 • t and r 3 = 7 • t x f. Each 
of these choices will generate different, nontrivial spin effects. Thus we believe that the 
essential feature of nontrivial spin transport should not be unique to our choice of boundary 
condition. Another direction of possible future work is the generalization of our approach 
to higher-dimensional spaces, as has been done with the nonrelativistic theory. 
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